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Abstract 

The elliptic flow v>2 of thermal photons at midrapidity in Au+Au collisions at -\Jsmn = 200 GeV 
is predicted, based on three-dimensional ideal hydrodynamics. Because of the interplay between 
the asymmetry and the strength of the transverse flow, the thermal photon V2 reaches a maximum 
at pi ~ 2GeV/c and the /^-integrated y 2 reaches a maximum at about 50% centrality. The pj- 
integrated V2 is very sensitive to the lower limit of the integral but not sensitive to the upper limit 
due to the rapid decrease in the spectrum of the transverse momentum. 



1. Introduction 

The deconfined and novel nuclear matter, the quark gluon plasma (QGP), has been expected 
to appear in relativistic heavy ion collisions. The observation of large elliptic flow of different 
hadronic species at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab, New 
York, is of special importance to confirm the the formation of QGR Unlike those bulk hadrons, 
photons are produced during the whole history of the evolution of the hot and dense matter. 
Moreover, the mean free path of photons is much larger than the transverse size of the bulk 
matter. So the produced photons pass through the surrounding matter without any interaction. 
As a result, thermal photons provide undistorted information on flow asymmetries not only from 
its surface but also from the inner of the hot, dense matter. 

Therefore the question such as what is the relation between the measurable elliptic flow V2 
of thermal photons and the evolution process of the expanding hot dense matter is of interest. 
This relation will serve as a direct bridge between the observables and the properties of the 
quark gluon plasma. A pioneering work has been done based on 2D ideal hydrodynamics, which 
has shown the transverse momentum and centrality dependence of the elliptic flow of thermal 
photons at midrapidity[l]. In this paper, the elliptic flow of thermal photons is calculated based 
on 3D ideal hydrodynamics. The paper is organized as following. In Sec. 2 we will briefly 
review the space-time evolution of the hot and dense matter using 3D ideal hydrodynamics and 
the basic formula for the production of thermal photons. In Sec. 3, we will show our results on 
the transverse momentum and centrality dependences of thermal photons V2 in Au+Au collisions 
at yfsNN = 200 GeV. Section 4 is devoted to discussion and summary of our results. 
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2. Bulk evolution, thermal photons, and elliptic flow 



In our calculation, a full 3D ideal hydrodynamic calculation is employed to describe the 
space-time evolution of the hot and dense matter created in Au+Au collisions at RHIC. The 
impact parameters corresponding to different centralities in Au+Au collisions at RHIC, which 
are estimated with a Glauber model, are 3.2, 5.5, 7.2, 8.5, 9.7, 10.8, and 11.7 fm for 0-10%, 
10-20%, ■ ■ • and 60-70% centrality, respectively. The local thermal equilibrium is assumed to be 
reached at the initial time To = 0.6 fm/c. The critical temperature of a first order phase transition 
between the QGP phase and the hadron phase is fixed at T c = 170 MeV. The initial flow is 
taken to be Bjorken scaling solution. So the transverse flow is vanishing at tq and dynamically 
generated by the pressure gradient. More details can be found in Q2|, l3[| . 

Transverse momentum spectra of thermal photons can be written as 
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with r(£"\r) being the Lorentz invariant thermal photons emission rate, d 4 x = Tdrdxdydr\ 
being the volume-element, the integration is done from the initial time to to the freeze-out, and 
E* = p^Ufj the photon energy in the local rest frame. More details can be found in 
The elliptic flow is quantified by the second harmonic coefficient V2 

f d(/>cos(2(l))d 3 N/dyd 2 p t 
Vl{P ^ y) = fdHWyd^ ' (2) 

where <p is the azimuthal angle of photon's momentum with respect to the reaction plane. 

The p, dependence of the triple differential spectra is strongly affected by the flow u through 
the argument E* = p^u^ in the photon emission rate. The azimuthal asymmetry of the transverse 
components of the flow obviously results in an anisotropic momentum distribution, which gives 
a finite V2- Therefore both strength and the anisotropy of transverse flow velocity are important 
to generate the elliptic flow of thermal photons. The two key features of QGP, mean radial flow 
(v r ) and the mean anisotropy of flow ^v^ )' aie defined as 
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<v r > = (V^)><r>(|^|), (3) 



where (■ ■ ■) stands for energy-density-weighted space-time average, v x and v v are the flow ve- 
locity components along jc-axis and y-axis, respectively. From Table Q] we see that the mean 
radial flow increases with centrality from 0-10% (b = 3.2 fm) to 20-30% (b = 7.2 fm), then 
decreases from 20-30% (b = 7.2 fm) to 60-70% (b — 1 1 .7 fm). On the other hand, the mean flow 
anisotropy increases with centrality from 0-10% to 60-70% monotonically. The increase of the 
average radial flow between 0-10% and 20-30% can be understood due to the fact that the initial 
shape of matter becomes more asymmetric, leading to flow asymmetry, as discussed earlier. The 
decrease beyond 20-30% is due to the decreasing life times, which prevents large flows from 
developing. 



3. Results 

In Fig. [TJ the V2 of thermal photons at midrapidity y = 0in0</?T<6 GeV/c is shown for 
various centralities from to 70% in Au+Au collisions at -\/Jnn = 200 GeV. The solid lines from 
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Table 1 : The strength and the anisotropy of transverse flow at each centrality. 

0-10 10-20 20-30 30-40 40-50 50-60 60-70 

3.2 5.5 7.2 8.5 9.7 10.8 11.7 

0.114 0.122 0.123 0.117 0.109 0.0959 0.0804 

0.0417 0.103 0.154 0.188 0.212 0.222 0.240 
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Figure 1: (Color online) Midrapidity v>2 for thermal photons in Au+Au collisions at ^snn 
various centralities from to 70% in < pj < 6 GeV/c. 



200 GeV is shown for 



bottom to top refer to centralities 0-10%, 10-20%, 20-30%, 30-40%, and 40-50%, respectively. 
The dashed lines from top to bottom at pj — 2 GeV/c refer to the centralities 50-60% and 60- 
70%. 

For each centrality, the thermal V2 increases then decreases with increasing p, and a peak 
appears at p t ~ 2 GeV/c. This p t dependence is consistent with the prediction based on 2+ ID 
hydrodynamics | lj and explained as the weak transverse flow at the early stage. 

The centrality dependence is shown clearly in Fig. [2] where the /^-integrated \'2 at y — 
is plotted as a function of impact parameter b. The pj -integrated v'2 reaches maximum at about 
50% centrality, due to the interplay between the asymmetry and the strength of the transverse 
flow, as shown in Tab. I. 

In Fig. [2] the dependence of /^-integrated range is presented with various types of curves. 
The /^-integrated V2 is not sensitive to the upper limit, as one can see that the curves from 
integrated range (0,3) and (0,6)GeV do not differ. But it is very sensitive to the lower limit, as 
shown in Fig. [2] where three lower limits of /^-integral, 0, 0.5 and lGeV/c, are used. The reason 
is that the pj spectrum of thermal photons decreases so rapidly|4j that the V2 weight at low pj 
spectrum plays an very important role. 



4. Conclusion 

Based on a three-dimensional ideal hydrodynamics, we calculated the elliptic flow of ther- 
mal photons V2 from Au+Au collisions at y/s^N = 200 GeV. Due to the interplay between the 
asymmetry and the strength of the transverse flow, thermal photon V2 reaches maximum at pi ~ 
2GeV/c and /^-integrated V2 reaches maximum at about 50% centrality. 
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Figure 2: (Color online) The px-integrated V2 at y = is plotted as a function of impact parameter b. The dependence of 
px-integrated range is presented with various types of curves. 



The insensitivity of the /^-integrated \>z to upper px-integral limit is very useful. In fact if we 
include all sources of direct photons, the px-integrated y 2 will not differ much, for two reasons: 

1) The yield of direct photons at higher pj is small because the pj spectrum of direct photons 
decreases rapidly though all sources of direct photons are counted[4]. 

2) The v>2 is also very small at higher pj for obvious reason that the main contribution at high p-f, 
the leading order contribution from nucleon-nucleon collisions, has a vanishing V2. 

So we can predict the /^-integrated V2 are dominated by thermal photons. 

The /?x-integrated V2 is very sensitive to the lower integral limit. Experimentally, the smallest 
measurable pj is up to the detector. Above this value, experimentalists can still vary the lower 
integral limit of pi to check how reliable the photon emission rate formula is. Theoretically this 
formula is very close to non-perturbative region and not reliable. 

Though 3D hydrodynamics provides similar results as from 2D hydrodynamics at midrapid- 
ity, the rapidity dependence will be certainly different, which will be presented later. 
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